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ABSTRACT 

The ultracompact binary systems V407 Vul (RX J1914.4+2456) and HM Cue 
(RX J0806. 3+1527) - a two-member subclass of the AM CVn stars - continue 
to pique interest because they defy unambiguous classification. Three proposed 
models remain viable at this time, but none of the three is significantly more 
compelling than the remaining two, and all three can satisfy the observational 
constraints if parameters in the mod els are tuned. One of t he three proposed 



models is the direct impact model of iMarsh fc Steeghd (120021 ). in which the ac 



cretion stream impacts the surface of a rapidly-rotating primary white dwarf 
directly but at a near-glancing angle. One requirement of this model is that the 
accretion stream have a high enough density to advect its specific kinetic energy 
below the photosphere for progressively more-thermalized emission downstream, 
a constraint that requires an accretion spot size ~1.2 x 10^ km^ or smaller. Hav- 
ing at hand a smoothed particle hydrodynamics code optimized for cataclysmic 
variable accretion disk simulations, it was relatively straightforward for us to 
adapt it to calculate the footprint of the accretion stream at the nominal radius 
of the primary white dwarf, and thus to test this constraint of the direct impact 
model. We find that the mass flux at the impact spot can be approximated by 
a bivariate Gaussian with standard deviation = 164 km in the orbital plane 
and a0 = 23 km in the perpendicular direction. The area of the the 2a ellipse 
into which ~86% of the mass flux occurs is r oughly 47,400 km^, or roughly half 



the size estimated by lMarsh &: Steeghd (120021 ). We discuss the necessary param- 



eters of a simple model of the luminosity distribution in the post-impact emission 
region. 



Subject headings: accretion — binaries: general — hydrodynamics — novae, 
cataclysmic variables — stars: individual(RX J1914. 4+2456) — white dwarfs — 
X-rays: binaries — X-rays: stars 
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Introduction 



Understanding tlie exact pliysical nature of tlie ultracompact binaries (UCBs) V407 Vul 
(RX J1914.4+2456) and HM Cnc (RX J0806.3+1527) has been the goal of extensive studies 
since these two closely related systems wer e discovered. The observational results tightly 
constrain the universe of possible solutions (ICropper et al.ll2004j : iBarros et al.l 120071 ): 



1. The "on/off" X-ray light curves folded on phase look nearly identical for the two 
systems, rising steeply to maximum, declining by half over the width of the pulse, and 
then dropping steeply again to zero 

2. The single observed and so presumed orbital periods are 321.5 s (5 .36 min) for HM 



Cnc flMotch &: Haberllll995l:ICroDper et al 



min) for V407 Vul jRamsay et al. 



2000 



1998 



Steeghs et all 120061 ). 



Ramsav et al.ll2002h and 569.4 s (9.49 



3. The observed per i ods are decreasing with time (jStrohmavei 



2002 



Ifiakala 



et al. 2003 



Strohmayen |2004J : iRamsay et al.l l2005l : iRamsay. Cropper, fc Hakala i2006l ) , consistent 



with the rate predicted for angular momentum losses by gravitational radiation, 

4. There is an apparent phase difference in the optical and X-ray light curves in both 
systems in which th e optical light curves lead the X-ray light curves by ~0.2 in phase 
tearros et al.ll2007h . 



5. No hard X-rays are detected. 



6. There is long-term variability in the mean luminosities of the systems ( llsrael et al. 



19991 : iRamsav et al.ll2002l : llsrael et al.ll2002[ ) 



7. Flickering, if present in the optical light curves, is of small amplitude, and suggests 
that th e optical light is no t dominated by radiation coming directly from the accreting 
region (IBarros et al.ll2007l ). 



Of the several models proposed to explain the early observations of these systems, three 
remain viable at this time: (i) the polar model, (ii) the direct impact (DI) model, and (iii) the 
unipolar inductor (UI) model. In all three of the currently viable models both the primary 
and secondary star are thought to be white dwarf stars, but only for the polar and DI models 
does the secondary fill it's Roche lobe. Each of the three can explain the observations with 
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varying degrees of success, but some fine tuning is required, and in the case of UI the physical 
process is unfamihar and unproven in a stellar context. 

Historically, the system V407 Vul was discovered first and suggested to be a dou- 
ble degenerate "soft" int ermediate pola r (IP) with a white dwarf spin period of 568 s 
(IMotch fc Haberll Il995l ) . ICropper et al.l (119981 ) noted that no other known IP displayed 
zero X-ray luminosity for ~l/2 the observed period, and they proposed that the system was 
inste ad a double degenerate analog of the pol ars. HM Cnc was identified in the following 
year (llsrael et al.lll999l : iBeuermann et al.lll999l ). and in 200 2 two separate groups suggested 



indep endently that this system is a near twin of V407 Vul (IRamsay et al.ll2002l : llsrael et al. 



20021 ) ■ In contrast, the binary ES Get with an orbital period of 10.3 min, only slightly longer 
than 9. 5 min for V407 Vul, shows clear signs of an ac cretion disk, and is a typical AM CVn 
binary (jWarner fc Woudtl l2002l : ICropper et al.ll2004l ). The IP model is unlikely to apply 
to these systems, in that the pulse shape of the X-ray light curve, the lack of hard X-rays 
and emission lines, and the infrared colors that se emingly exlude a mai n-sequence secondary 
are all difficult to explain within the IP model (ICropper et al.ll2004j). We do not discuss 
this ra odel further, but the interested reader may wish to consult iNorton. Haswell. &: Wynn 
J2004j ). 



The polar model for the V407 Vul stars proposes that a Roche-lobe filling secondary 
transfers He-rich plasma to the magnetized and synchronized primary star, channeled by 
the magnetic field onto the nearest magnetic pole which is fixed in the co-rotating frame. 
Orbital motions take the accretion spot and column in and out of view, producing the X-ray 
lightcurves. The phase-offset optical variations are thought to arise from reprocessing on the 
face of the secondary star. Problems with this model include that emission lines and polar- 
ization are common in polars but not observed here, and mass transfer from a lower-mass 
secondary to a higher-mass primary sho uld increase the orbital period. In regards t o this 
last point, however, two separate groups (jPeloye fc Taanj|2006l : iD'Antona et al.l 120061 ) have 
shown that secondaries that are still contracting and non-degenerate can yield decreasing 
orbital periods while the mass transfer is in early phases, and can exist in this state for 10^ 
to 10^ years. This is sufficient t ime for the requireme nts of observed space densities based 
on population synthesis models (INelemans et al.ll200ll ). 



The DI model proposes that the two stars are clos e enough that the accr etion stream 
impacts the surface of the primary white dwarf directly (IMarsh fc Steeghsll2002l ) because the 
stellar radius is lar ger than the "per i astron " of a ballistic stream originating at the inner 
Lagrange point LI. iMarsh fc Steeghd (120021 ) argue that in the DI model the accretion stream 
will be sufficiently narrow (width ~ 1Q~'^Rq) that the density of accreted material will be 
high enough to penetrate below the photosphere and thermalize to soft X-rays, similar to 
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the blob accretio n model proposed for some polars (iKuijpers fc Pringld Il982l : iFrank et al. 
19881 : iKina I2OOOI ) ■ For this to apply in the DI model, the surface area of the accretion spot 
must be a fraction / ~ 10~^ of the surface area of the star or smaller, which translates to an 
area smaller than ~ 1.2 x 10^ km^. The initial motivation for this project came from Pasi 
Hakala, who suggested to one of us (MW) that it would be useful to use our SPH code to 
determine accurately the accretion spot size. If a hydrodynamic sir nulation were to find that 
the footprint of the accretion stream is substantially larger than the lMarsh &: Steeghd (120021 ) 
estimate, then that would be a serious problem for the DI model. Of course, it is a one-sided 
test, and passing the test does not mean that the DI model is correct, just consistent with 
the physical constraints of the blob accretion model. 



The unipolar-inductor (UI) model proposed by IWu et al.l (|2002|) a n d discussed most 



recently in th e context of V407 Vul and HM Cnc by iDalFOsso et all (120061 . 120071 ) and 



Ramsay et al.l (120071 ) is the stellar analog of behavior observed in the Jupiter-Io system. 
As lo moves through Jupiter's magnetic field, charge separation occurs, and because the 
environment is a weak plasma, an emf is established from the far point of lo as observed 
from Jupiter, along the magnetic field lines in both directions t o Jupiter's atmosphere in 



both hemispheres, and from the planet back to the near side of lo (IGoldreich fc Lynden-Bell 



19691 ). Resistive heating in the region between the two spots in each pair resulting from the 



(Clarke et al. 


19961 


Wu et al. 


(2002) 



could explain the observations - certainly a fascinating and clever idea, even if eventually 
determined not to be applicable in these stars. 

In a stellar context, the rate of energy dissipation in the UI m odel can be of or der Lq 
or higher, and with the predicted hot spot area ~ 8 x 10^ km^, the IWu et al.l (12002! ) model 
predicts a 50 eV temperature co nsistent with the 55 eV observed from ROSAT and ASCA 
spectral fits ( Ramsay et al.|[200ol ). The hotspots are on the plane that passes through both 
the line of centers and the orbital rotation axis, and the source of the optical light modula- 
tion is again assumed to be reprocessed radiation from the heated face of the secondary. An 
important point in this model is that the secondary is not Roche-lobe filling, and there is no 
significant mass transfer. To first order the timescale for orbital evolution is that which grav- 
itational radiation predicts, consistent with the observed rates of period change (IStrohmavei 



2002 



Hakala et al. 



2003 



Strohmaverll2004 iRamsay et al.ll2005t iRamsay. Cropper. &: Hakala 



20061 ). iRamsay et al.l (120071 ) note that the UI model predicts radio emission, which they 
detected at the 5.8-0" level along the line of sight to HM Cnc, but which was not detected in 
observations of V407 Vul. The major problem with this model is the phasing of the optical 
versus X-ray light curves, because the X-ray spots are on the meridian facing the secondary, 
this model predicts a phase offset of ~0.5, whereas the observed phase offset is ~0.2 for both 
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V407 Vul and HM Cnc flBarros et al.l 120071 ). In addition, Marsh fc NelemansI feoosh sue.- 
gest there are serious problems with synchronization timescales, but iDall'Osso et aLl (120071 ) 
argue that Marsh & Nelemans result rests on an estimate of the luminsoty of V407 Vul 
which is too high by a fact or of ~ 10^. However, for the UI model to explain both sys- 



tems, IDairOsso et al.l (120071 ) find that the primary white dwarf in HM Cnc rotates nearly 
synchronously but slightly slower than the orbital period, whereas the primary in V407 Vul 
must have a rotation frequency which is roughly 10% higher than the orbital frequency. The 
authors do not suggest how prior binary evolution without mass transfer can yield such a 
rapidly-rotating primary. 

In all three of the models above, the on/off behavior of the X-ray light curve has been 
i nterpreted to mean that t he X-rays must be emitted from a geometrically compact region 
(jCropper et al.l 119981 . 120041 ). For the UI model, this region is the footprint of the magnetic 
field l ines which ca r ry cu rrent from the secondary (non- magnetic) to primary (magnetic) 
star. 



Barros et al.l (120051 ) analyzed this model and were able to place strict constraints 



on the allowable parameter space based on geometrical arguments involving the size of this 
footprint. Here, we present the results of a hydrodynamic simulation of V407 Vul to estimate 
the area of the footprint of the accretion stream, to determine if the calculated spot size is 
small enough to satisfy the DI model requirement. In §2 we discuss the changes made to 
our SPH code FITDlSK to pursue this project and the system parameters assumed. In §3 
we discuss the characteristics of the simulated accretion stream, and in §4 we present our 
primary results. In §4 we discuss our results in the context of the existing literature on V407 
Vul and HM Cnc, and in §6 we present our conclusions. 



Numerics 



2.1. The Code 



We use the method of smoothed particle hydrodynamics, a Lagrangian-based technique 
that is well-suited fo r astrophys i cal problems with unu sual geometries and vacuum outer 
boundary conditions (lLucylll977l : lMonaghanlll992l . 120051). Our cod e FITDlSK was deyolped 



to model accretion disks in cataclysmic variables (ISimpsonI Il995l: ISimpson &: WoodI Il998 
Wood et allbood : IWood. Dolence fc Simpsonlbood : IWood fc Burkell2007h . and the details of 



the numerics can be found in those publications. In brief, FITDlSK solves in the inertial 
frame the equations of compressible, viscous hydrodynamics subject to the external gravita- 
tional field of two point-mass stars. Particles all have the same effective radius (smoothing 
length h) but can have different time step sizes as dictated by local conditions. In our typ- 
ical simulation runs, we build up a disk by injecting ~2000 particles per orbit at the LI 
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point until we reach the target number of particles in a particular run (~10^). A particle is 
considered accreted onto Mi when it is closer than the specified radius of the star. When 
this occurs (or accretion onto M2 or ejection from the system), a replacement particle is 
immediately injected at the LI point. Thus, the simulation in time settles into a dynamical 
equilibrium state. 



Our goal with the current project is to test the requirement of the DI model of lMarsh fc Steeghs 



(120021 ) that the accretion spot on the surface of Mi has an area A < 1.2 X 10^ km^ such that 
the stream density is sufficient to advect and release the stream kinetic energy below the 
photosphere, which is required to avoid hard X-ray emission. For this preliminary test, then, 
our main requirement is that the accretion stream behave as a fiuid from the LI region to the 
surface of the primary. We keep track of where the trajectories of the simulation particles 
intersect the mathematical surface of the primary, and use the resulting integrated mass-fiux 
distribution to estimate the effective area of the acccretion spot. Our goal here is not to 
model in detail the hydrodynamics of the impact region - we are only beginning to address 
that far more ambitious calculation, and it will be the topic of a future publication. Our 
current results can be considered to provide the minimum estimated effective area for the 
accretion spot, as all physical effects we've neglected could only act to increase the effective 
area. 

To explore this problem, we made two substantial changes to our base code. First, 
we converted the code to make the calculations in the co-rotating frame instead of the 
inertial frame used in FITDlSK up to this point, and second we implemented a new particle 
injection scheme which gives a far more physical model of the fiow dynamics through the 
inner Lagrange point LI. 

Our original code calculates in the inertial frame, which generally simplifies the analysis 
and visualization steps, but for this project the analysis is simplified if the calculations are 
completed in the co-rotating frame with the implementation of coriolis forces. We establish a 
Cartesian system, with the center of mass at the origin. The x-axis is along the line of centers 
of the two stars, with the donor star orbiting in the +y direction and the orbital angular 
momentum vector in the +z direction. Within this system, the stars, stream, and impact 
spot are stationary features. Therefore, the approximate point of impact for a given particle 
can be easily calculated by linear interpolation along its trajectory between its positions 
at time steps and t" bracketing the primary surface radius. Also, in a steady state 
it becomes meaningful to examine the cumulative distribution of the particle mass fiux as 
a function of position on the surface of the primary as a means to estimate the effective 
accretion spot size. 

The mass transfer stream is created ah initio by a Monte- Carlo-based particle injection 
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technique. Particles are injected at random locations in a small region near the LI point. 
In the original code, the region is a box located on the primary side of the LI point, and 
the LI point is at the center of one face. Particles are injected with velocity vectors cal- 
cu lated assuming soun d-speed flow through LI, and we naturally recover the trajectories 



of iLubow fc Shul (I1975I ). In order to resolve the accretion spot, we used an effective scaled 
particle size of 10 km, or /i ~ 10~^ in simulation units where the stellar separation a = 1. 
With particles this small, we found the particle trajectories were ballistic and not behaving 
as a fluid (cf. Cash 2002, 2004). The solution we arrived at was to move the injection region 
to the secondary side of LI, and to implement a reflecting boundary plane approximately 
30/i from the LI point. We inject particles with zero velocity in this region, and hence create 
a crude pressure-supported atmosphere that naturally leads to a physical flow distribution 
through the LI region into the Roche potential of the primary. To avoid spuriously generated 
sound waves when injecting particles, we find that preferentially injecting particles near the 
base of the atmosphere, where the density is significantly higher, is an effective technique 
since viscous interactions tend to quickly damp any excited modes. This technique provides 
the most accurate flow dynamics through the LI region obtainable within the SPH approach. 

The Reynolds number is deflned as 

v.L 



Re 



5 



where Vs is the mean fluid velocity L is the characteristic size of the region, and v is the 
dynamic fluid viscosity. We can es timate the Reyn olds number in the LI region in our 



simulation by using the results of (iMonaghad 119971 ) who flnds that the kinematic shear 
viscosity within the SPH formalism can be approximated as 

15 

u = —Kvsigh (2) 

where is a constant {K ~ 0.5) and v^ig = Ci + Cj is the signal velocity between particles 
i and j. As noted above, we place our reflecting surface 30h behind the mathematical LI 
point, which is roughly the width of the stream at the LI point (see Figure 1), and we use 
this as the characteristic size L (in pipe flow, the width of the pipe is taken as L). The mean 
fluid velocity is f s ~ at the reflecting wall, and Vg ~ as the flow passes through the LI 
point, so we use Vg = Cs/2 in the estimate of the Reynolds number. Then from equations (1) 
and (2) and using i^' ~ 1/2, t>sig = q + Cj ~ 2cs, and L ~ 30h, we have 

v.L (c,/2)(30/i) 

-n^e = ~ -j-E-^ — -— — ~ 112. (6) 

V ^(l/2)(2c,)/. ^ ^ 

This value of the Reynolds number is roughly 3-4 orders of magnitude lower than the typical 
transition to turbulent flow, so we expect the fluid near the LI region to be non-turbulent. 
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Furthermore, we note that iDeloye et aLl (120071 ) in their recent study of the structure and 



evolution of the mass donating secondary stars found that the atmospheric structure of 
these stars just after mass transfer is initiated is characterized by radiative stabihty at the 
photosphere, a thin and weak sub-photospheric convection zone, and then another radiative 
zone that extends to the core. 

As discussed above, particles that are accreted onto the primary are immediately rein- 
jected in the LI region so that a constant number of particles is maintained in the simulation, 
and a dynamical equilibrium flow is established. For the simulation presented here, a steady- 
state mass transfer rate is established less than one orbit after the number of particles reaches 
its preset maximum value. 



2.2. The Model 



For the current study, we adopt the system parameters suggested by iMarsh &: Steeghs 



(120021) for V407 Vul: Porb = 9.5 min. Mi = O.SMq, Ms = O.IM©, and Ri = O.O144i?0. The 
orbital separation is calculated from Kepler's third law and is taken to be constant. Distances 
are scaled so that the orbital separation is unity in dimensionless system units. Because the 
mass of the SPH particles formally drops out of the equations (see SW), no particular mass 
transfer rate is assumed a priori. The total number of particles is set to 86,400. Particles are 
given an initial tempera ture of 20,000 K, appropriate for t he heated face o f the d onor star 



(IMarsh fc Steeghsll2002l ). The viscosity coefficients for the iLattanzio et al.l (119861 ) artificial 
viscosity prescription are a = 1.0 and (3 = 0.5, although viscosity is expected to play a 
negligible role in the dynamics of a stream of plasma in free-fall. An ideal equation of state 
is used, with 7 = 1.01. The smoothing length was set to a value expected to be small with 
respect to the impact spot, which in our case was h = 1.1493 x 10""^ in simulation units, or 
h ^ 10 km. Tests were done with slightly larger h and the result did not change significantly. 
With a smoothing length of this size, we were forced to use 6,000 time steps per orbit with 
three levels of refinement (i.e., 6tQ = Prmorfe/6000, 6ti = 51^/2, or 5t2 = Sto/A) as described 
more fully in SW. We note that the time-of-fiight t/ from LI to impact can be estimated 
using Kepler's third law with a' ~ a/2, and then the time of fiight from LI to impact is 
roughly tf ~ (l/2)(a'/a)^/^Porb ~ Porb/5 ~ 1-2 min, where the factor of 1/2 is the result of 
the particle completing only half the orbit. 

While we make an attempt to accurately calculate the stream profile from the LI region 
to impact spot, we make several simplifying assumptions. First, we assume that rotation of 
the secondary and surface fiow on the secondary star were negligible. This is likely to be a 
reasonable assumption as the secondary is almost certainly co-rotating at the binary period. 
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We further assume for the purpose of these calculations that in the absence of magnetic 
fields the accretion spot size will not be greatly enhanced by interactions with the fluid 
already at the impact point, and we assume that each particle impacts at the point where 
it's trajectory takes it through the radius of the primary. This assumption should hold to 
a good approximation, because the fluid in the accretion stream is highly supersonic when 
it impacts the surface, and in the absence of significant surface magnetic fields, the surface 
fluid is unable to apply any force on the incoming fluid before impact with the exception 
of (negligible) radiation pressure. Finally, in all calculations, radiative effects and magnetic 
forces are explicitly ignored. 



2.3. Analysis 

As described above, our technique lends itself to examining the cumulative distribution 
of particles over a long span of simulation time. Once our simulations have reached dynamical 
equilibrium, the probability distribution function of particles is constant with time. By 
utilizing a Monte-Carlo approach to particle injection, we ensure that the particular positions 
of the particles changes over the course of the simulation. We further ensure a random 
selection of particles by sampling the simulation every GO*'^ time step so that a particular 
particle does not contribute to the number density of particles in any location more than once. 
Note that particles within the Roche-lobe of the secondary could potentially be sampled 
many times by this technique, but this is of no concern here. With these aspects in mind, a 
two-dimensional Eulerian mesh (with small thickness) is created that slices through a region 
of interest and the number of particles that lie within each grid zone is summed over all 
of the time steps included by the above criterion. In this way, we effectively increase the 
resolution of our results and reduce the noise inherent to the SPH technique (see also Wood 
et al. 2005). 

To characterize the impact spot, many of the above ideas apply. The key difference is 
that we record the impact position of every particle accreted, as opposed to sampling every 
60*'* time step. This is because each particle can only cross the surface of the primary one 
time by design, and so the Monte-Carlo procedure for particle injection, coupled with the 

significant mixing within the Roche-lobe of the secondary, serve to randomize the particle 
impact points. Also, since the impact spot lies on the curved surface of the primary, we lay 
our mesh on this curved surface and utilize the angular position of the impacts to bin them 
into grid zones. While this creates zones of different linear dimensions, the effect is negligible 
since the impact spot lies in a very small region centered on the orbital plane. Once the 
particles have been binned, the relative positions of the grid zones are converted back to 
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linear dimensions so that the distances presented below are actually arc lengths. With this 
complete, the distribution is fit with a bivariate Gaussian using a nonlinear least squares 
analysis. 



3. Stream Characteristics 

In this section, we briefly recount the results of several previous studies which explored 
the dynamics of the mass transfer stream in accreting binary stars. We focus on two regions, 
namely the LI region and the impact spot on the primary, the former primarily for validation 
of our technique. 



3.1. The LI Region 



In the LI region. iLubow &: Shij (Il975l . Il976l ) showed via semianalytic methods that the 
density cross section of stream should Gaussian, and decrease in overall density as where 
r is the distance from the LI point along the center of the stream. They also showed that 
the stream is deflected by an angle 6 ~ 20.8° from the line connecting the centers of the 
stars (for mass ratio q = 0.2). Later numerical work by (jOka et al.ll2002l ). which accounted 
for the surface flow on the secondary, found that the stream was approximately Gaussian 
about its center (slightly bow shaped convex towards the trailing side), decreased as r~°-^, 
and had a deflection angle smaller by roughly a factor of two compared to Lubow & Shu. 



3.2. The Accretion Stream Impact Region 



As the stream material accelerates toward the the primary, it is stretched by tidal 
for ces which act t o nar row the stream, especially in the z-direction. This was predicted 
by iLubow fc Shul (Il976l ). who showed that the vertical stream width at closest approach to 
the primary should be of order the hydrostatic scale height of a disk at that same radius. 
The pressure scale height of a hydrostatic disk at r = i?i from the center of the primary is 
approximately given by 

/ kTR^, \ , 
— — ^ ~ 6 X 10"^ 
\GMimp) 



in simulation units. iLubow fc Shul (119751 ) find that the stream width in the orbital plane 
should be roughly a factor of five larger than the vertical width at closest approach, or 
~ 3 X 10~^ in simulation units. While in our case the stream is intercepted by the surface of 



the primary before it reaches its closest approach, we take the se values as approximate, a nd 
indeed find our results to be in good agreement with those of iLubow fc Shul (Il975l . Il976l ). 



Results 



In Figure 1, we show the trajectory of the stream by plotting the stars and the SPH 
particles at a single time step viewed at orientations perpendicular and parallel to the orbital 
plane. As previously noted, t he impact spot on the p rimary white dwarf is on the far side 
as viewed from the secondary (IMarsh &: Steeghsll2002l ). 



As a proof of concept, we begin by analyzing our results for the fiow in the LI region. 
Our results are shown in Figure 2 as density contours of a slice centered on the orbital 
plane. We find that the stream is indeed approximately Gaussian distributed about its 
center line, decreases in density as approximately r~°'^ , and i s defiected by an angle 9 ~ 21°, 



in reasonable agreement with the iLubow fc Shul (1 19751 . Il976l ) results. 



In Figure 3, we show contours of constant mass fiux through the surface of the primary 
and the best fit obtained via nonlinear least squares analysis using a bivariate Gaussian 
function. For the purposes of fitting, a coordinate system with origin at the center of mass 
fiux is employed. The mass fiux, $, as found by our model is of the form 



se) = N^exp 



2 \ a! + a2 



(4) 



where iV$ is a constant scaling factor, = (pRi is the arc length in the azimuthal direction, 
S0 = 6Ri is the arc length in the polar direction, and and ao are the standard deviations 
in these directions. In these coordinates, increases in the —x direction and se increases 
in the +z direction. The result shown in Figure 3 is normalized to unity at the origin. The 
best fit model gives = 0.0019 and ag = 0.00026 in simulation units, or cr^ = 164 km and 
ag = 23 km in physical units. The eccentricity of the resulting ellipse is e = ^1 — {ag/a^Y ~ 
0.99. 



5. Discussion 

Looking again at Figure 1, it is apparent that the supersonic accretion stream impacts 
the white dwarf surface with a nearly horizontal trajectory. The altitude angle of the stream 
as viewed from the impact point is <20°). A very approximate mental picture of the situation 
can be achieved by considering a narrow, high- velocity jet of water impacting the surface of a 



- 12 - 



pool with a similarly shallow altitude angle, and in fact the experiment is so easily available 
we recommend the interested reader actually try it. In the stellar case, the shock-heated 
surface fluid will be promptly carried away from the impact point, and the accretion stream 
will always be impacting fresh material that is entrained from the surrounding T ~ T^s 
atmosphere. There will be strong radial flows towards the core of the impact spot which will 
serve to advect on a somewhat slower time scale the speciflc kinetic energy of the accreting 
material in the lower-density outer volume of the stream. 

We characterize the shape of the impact region as an 'W-ci" ellipse with semi-major and 
semi-minor axes given by for example a = and b = (Xe for the 1-a case, and where N need 
not be an integer. If we now assume our coordinate system is centered on the ccntroid of the 
mass-flux distribution, and that the distributions along both the and 9 axes are normal 
gaussians, then the constraining equation for the bivariate normal density distribution for 
the N-a contour is (for example, see Kawka et al. 2004, §4.2) 

or 

\ 2 / \ 2 



+ - <Ar^. (6) 
(TcpJ \crej 

The joint probabihty of finding an object within the N-a ellipse is then the integral of the 
probability density function inside the ellipse: 



1 rN I'^/IS'—U-' 

P(Na)^— e-"'/M e-'^'/^dvdu, (7) 

27r J_N J-^N'^-ui 

where u = s^/a^ and v = s^/ag. This yields for the 1-, 2-, and 3-a countours: 

P{la) = 0.36 (8) 
P(2a) = 0.86 (9) 
P(3a) = 0.99. (10) 

Thus approximately 86% of the mass in the accretion stream impacts a spot of area — 
47, 400 km^, or roughly a fraction /a^ = 3.8 x 10"^ The S-a profile (99% of the mass fiux) 
yields an area A^^ — 107, 000 km^, or roughly a fraction f^a — 8.5 x 10~^, or about the size of 
Iceland on a white dwarf of radius ~1.5i?0. Our numerical results confirm that the footprint 
of the accretion stream is remarkably compact, and is consistent with the requirement of the 
DI model we set out to test. 



We can determine whether the accretion stream penetrates below the photosphere by 
comparing the ram pressure of the accretion stream to the gas pressure in the white dwarf 
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atmosphere. If we take a conservative estimate for the accretion rate M = 10~^° M0 yr~^ 
(jPeloye et aLl 120071 ). then we find that the distribution of the stream ram pressure in the 
radial direction has the same form as the mass flux given above, with a peak value of 
3.4 X 10^ dynes cm~^, as shown in Figure 4 . White dwarf model atmosphere calculations 
JVennes fc FontainelEog^ : iKawka et allboOTh for a T^s ^ 20, 000 K, log^ = 7.75 white dwarf 
yield gas pressures at optical depths r = 1 and 10 of Pgas ~ 3 x 10^ and 2 x 10^ dynes cm^^, 
respectively. Thus, Pram ^ -Pgas, and more than 99% of the mass flux should advect below 
the photosphere where the bulk of its energy can be thermalized as in the blob accret ion 
model proposed for some polars (iKuijpers fc Pringldll982l : iFrank et al.lll988l : lKingll2000l ). 



Marsh &: Steeghd (120021 ) suggested that because the specific kinetic energy of the accre- 



tion stream enters at such a shallow angle to the surface, there should exist a high-luminosity 
region just downstream of the accretion impact point that emits in soft X-rays, and further 
downstream in the prograde equitorial flow one would expect th e spectral energy distr ibution 
to soften further towards mean stellar effective temperature. iRamsay et al.l (120051 ) report 
spectral softening as the flux declines to wards the off-phase i n V407 Vul. The mean effective 
temperature of HM Cnc is estimated by iBarros et al.l (120071 ) to be in the range 18,500 K to 
32,400 K (and probably similar for V407 Vul but currently unmeasurable because of the G 
star which is conincident with that binary). We assume T^s ~ 25, 000 K for the purpose of 
some simple estimates discussed below. The peak optical radia tion will thus lead the X-ray 
peak in orbital phase. Recent observations (IBarros et al.l 120071 ) indicate that both systems 
HM Cnc and V407 Vul are characterised by the opt ical pulses leading the X -ray pulses by 
about 0.2 in orbital phase, and that earlier claims (jlsrael et al.l l2003l . |2004| ) that the opti- 
cal and X-ray pulses were antiphased were the result of incorrect timing corrections for the 
optical data. 

We can crudely model the surface luminosity distribution using 4 zones: the first three 
are the accretion spot, the soft X-ray spot, the optical spot - each of which is assumed to be 
uniformly illuminated - and the forth is the remaining stellar photosphere. All three of the 
spots will be highly-elongated features, and very roughly we approximate them as ellipses 
with a common focu= (the center ot the impact .pot) and eccentridtja (see Figure 5). As the 
newly-accreted and shock-heated gas flows along the equator, shear mixing and high levels 
of turbulence will effectively spread the width of the accretion luminosity flow in latitude as 
the material propogates in longitude away from the impact spot. 



^We note that this approximation is simply the assumption that the area of the spots are proportional 
to the square of the length - the results would be essentially the same for scaled rectangular regions where 
all 3 aligned on one edge. 
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Because the X-ray flux is observed for 60% of the orbital period, the soft X-ray emitting 
region must have a linear extent ~10% the circumference of the star, or equivalently 0.1 in 
orbital phase. This corresponds to a physical length of s ~ 0.1(27ri?i) = 6100 km, which 
is roughly 9 times longer than the full length of the 2-cr accretion spot discussed above. 
Measuring from center of the accretion spot, the distance to the center of the X-ray spot will 
be roughly 0.05 in phase, or 3050 km. Because for high eccentricity ellipses the distance from 
the focus to the center of the ellipse is ea ~ a, for a uniformly-illuminated optical spot to 
lead the X-ray spot by 0.20 in phase its center would be at a phase angle of 0.25 as measured 
from the accretion spot, which then implies that the optical spot extends for ~0.5 in phase 
around the primary white dwarf (Figure 5). In the limit where the area scales as the square 
of the ellipse length minus the area of the next smaller ellipse, the ratio of areas of the 3 
zones is then 1:80:1920, and the ratio of the areas of the X-ray to optical spots is 1:24. In 
the limit where the areas of the zones scale simply as the length of the zones (constant width 
of order the impact spot width), the ratio of areas for the 3 zones would be approximately 
1:9:45. We note that the real emission regions will have monotonically-decreasing radiation 
fluxes as a function of distance from the impact spot, which then means that the above 
estimates are lower limits to the linear extents of the X-ray and optical zones. 



Barros et al.l (120071 ) reports that the peak-to-peak amplitude for HM Cnc is roughly 
25% in the optical. For a 2-a spot viewed in the center of the stellar disk, this implies a spot 
temperature which is some 16% higher than that of the surrounding photosphere (~29,000 
K). At an inclination oi i = 45°, the spot temperature needs to be some 40% higher than 
that of the surrounding photosphere (~35,000 K). The amplitudes in V407 Vul are diluted 
by the the G-star which is along the same line of sight and thus not known unambiguously. 

The DI model appears to be able to accou nt for the observ a tions of bo th systems V407 



Vul a nd HM Cnc, if we accept the results of iDeloye fc Taanj (120061 ) and iD'Antona et al. 



( I2OO6I ) that a contracting lower-mass secondary can yield a negative rate of period change 
and a steady accretion rate for a sufficient length of t ime that the obiects' space density 
is high enough to be consistent with the observations. ICropper et al.l (120041 ) suggests that 
the character of the X-ray light curves of these two systems is difficult to understand in the 
context of the DI model, where they propose that within the DI model one would expect an 
X-ray ligh t curve that showed a slow rise to maximum and a steep decline, opposite to that 
observed. iRamsay et al.l (120051 ) find that the spectral softness varies over the orbital phase, 
and interpret this to mean that the X-ray emitting region has some temperature structure. 
More work - both observational and theoretical - will be required to understand the nature 
of the luminosity distribution on the surface of the white dwarf primary within the DI model, 
and whether there is a unique solution that explains the set of observations. 
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6. Conclusions 



We have shown that for system parameters appropriate to V407 Vul, the impact spot 
size estimated with a hydrody namics simulation i s sma ll enough to be consistent with the 
requirement of the DI model of iMarsh &: Steeghd ((2002) ■ In summary, we find: 



1. The impact spot calculated for the system parameters of V407 Vul is very compact, 
with an effective area only 47,400 km^ or a fraction / = 3.8 x 10^^ that of the of the 
primary star assuming the 2-0" contour through which 86% of the mass flows. The 
3-0" profile (99% of the mass flux) encloses an area A^^r = 107,000 km^, or roughly 
a fraction /3a- = 8.5 x 10~^ of the stellar surface area. The mass flux at the spot 
is characterized by a bivariate Gaussian distribution with standard deviations a,/, = 
164 km and ag = 23 km. 

2. Within the DI model, the accretion stream impacts the surface nearly horizontally 
(altitude ^20°) and supersonically. The specific kinetic energy is advected and ther- 
malized below the photosphere to be reemitted downstream from the impact point in 
a very elongated region along the equator and this equitorial flow is characterized by 
decreasing localized temperature and therefore decreasing spectral energy distribution. 

3. The X-ray emitting region must extend for ~10% of the equitorial circumference to 
account for the visibility of the X-ray source for ~60% of the orbit. This region 
then should be centered roughly on phase 0.05, measured along the equator from the 
impact point. For the optical phase to lead the X-ray phase by 0.2, the optical region 
is centered on phase ~0.25, which then implies the region has a full extent of ~0.5 in 
phase, or half or more the circumference of the accreting white dwarf. 

4. The fluid temperature must decrease from ~50 eV (~500,000 K) to ~25,000 K in the 
time it takes the equitorial flow to bring the newly-accreted fluid roughly half-way 
around the primary. For a flow speed equal to the breakup velocity, this gives a time 
of 12 s, and an exponential cooling timescale of ~4 s. For a more reasonable flow speed 
that is one-tenth the breakup velocity, the cooling timescale is ~40s. Note that in this 
context "cooling" includes both radiative losses and entrainment of cooler gas from the 
surrounding photosphere. 



Of the three viable candidate models for the observed systems V407 Vul and HM Cnc, 
the DI model appears to require the least number of "just so" parameter choices, and to 
fit most easily the observational constraints. The major issue with the DI model is the 
observed negative rate of period change, but as discussed above finite-temperature models 



- 16 - 



of the secondary stars appropriate to these systems indicate overall contra ction which can 
yield stable mass transfer an d a decaying orbit for a sufficient span of time (jPeloye fc Taam 
20061 : iD'Antona et al.ll2006l ). The polar model is then also not excluded by the observed 
dP/dt rates. Within the polar model the accretion impact spot is likely to be located 
within 7r/2 radians from the stellar line of centers, and this makes the observed phasing 
problematic. The accretion stream trajectory in polars should be nearly vertical at impact, 
and because it is magnetically confined there should be a standing shock. The absence 
of a hard component in the spectrum as well as the absen ce of polarization are further 
long-standing problems with this model (ICropper et al.ll2004l ). The unipolar inductor model 
appears to allow consistency with all of the observations - including the reported detection 
of radio emission (IRamsay et al.ll2007l ) - except the relative phasing of the optical and X-ray 
light curves. In the UI model the ~50 eV hot spots on the surface of the primary white 
dwarf are found along the meridian that passes through the line of centers, and the optical 
light curve results from the heated face of the secondary star. This model nicely fit the 
old, incorrectly time-calibrated observations that the optical and X-ray light curves were 
anti-phased, but there appears to be no way for this model to explain the now currently 
accepted phase offset of 0.2 determined by iBarros et al.l (120071 ). 



In future work, we plan to model the hydrodynamics of the impact region, using an 
implicit technique for estimating the radiation field and energy exchanges between SPH 
particles in the fiux-limited diffusion approximation as described in lWhitehouse et al.l (120051 ). 
and to use simple ray-trace concepts to estimate the emergent radiation as a function of 
inclination angle and orbital phase. 



We thank Pasi Hakala for suggesting this project and Adela Kawka and Tom Marsh for 
useful discussions. We also thank the referee for suggestions which improved the presentation 
of the results. This work was supported in part by grant AST-0205902 to The Florida 
Institute of Technology. 
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Fig. 1. — Stream profiles in the planes parallel and perpendicular to the orbital (x-y) plane. 
Also, the horizon of the impact spot is shown as the dotted line. Note that the altitude of 
the incoming accretion stream at the white dwarf surface is <20° above the local horizen. 
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Fig. 2. — Equi density contours in the LI region. The density is normahzed to unity on the 
contour that crosses over the LI point. Note the different scales on each axis. 
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Fig. 3. — Contours of constant mass flux at the impact spot. The origin is set to the center 
of the fitted distribution (dashed hnes). The values shown have been normalized to unity at 
the central maximum value. Note the different scales on each axis. 
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Fig. 4. — Contours of constant ram pressure in the radial direction at the impact spot. The 
values shown are in units of 10^ dynes cm~^. 
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Fig. 5. — Schematic diagram of a crude 3-zone model. Assuming the impact spot, X-ray 
spot, and optical spot can all be characterized as uniformly-illuminated ellipses with the 
same eccentricity (e = 0.99 is shown), and using the observational constraints, wc find that 
the X-ray spot extends roughly 0.1 in phase, while the optical spot would extend ~0.5 in 
phase, or half the circumference of the primary white dwarf. 



